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We compute the chiral corrections to octet baryon axial currents through 0{p^) in heavy baryon 
chiral perturbation theory, including both octet and decuplet baryon intermediate states. We in- 
clude the latter in a consistent way by using the small scale expansion. We find that, in contrast 
to the situation at 0{p^), there exist no cancellations between octet and decuplet contributions at 
0(j)^). Consequently, the 0{p^) corrections spoil the expected scaling behavior of the chiral expan- 
sion. We discuss this result in terms of the 1/Nc expansion. We also consider the implications for 
determination of the strange quark contribution to the nucleon spin from polarized deep inelastic 
scattering data. 



PACS Indices; 12.39. Fe, 13.30.Ce, 14.20.Jn 
I. INTRODUCTION 

The chiral expansion of the octet baryon axial current 
J^5 has been a topic of ongoing theoretical interest for 
some time. At 0{p^), this current is parameterized by 
the well-known SU(3) reduced matrix elements D and F. 
The leading chiral corrections, which arise at 0{p^) con- 
tain chiral logarithms, which were first computed in Ref. 
10,^. Subsequently, the wave function renormalization 
correction was added in the framework of heavy baryon 
chiral perturbation theory (HBCPT) |^-|^, which pro- 
vides for a consistent power counting. While these cor- 
rections are large when only octet baryon intermediate 
states are kept inclusion of decuplet contributions 
produces sizeable cancellations, leading to a significantly 
smaller 0{p'^) effect Q. The origin of these cancellations 
may be explained by considering the large Nc expansion 
1^, as noted in the work of Refs. |0-|ll|. In terms of 
this counting, the 0{p'^) contribution are of order iVc, 
while the 0{p'^) loop corrections are nominally 0{N^). 
As shown in Refs. |p|-|l2| , however, a spin-fiavor symme- 
try arises at this order whose effect is to render the 0{p'^) 
loop effects of relative order N^- Thus, inclusion of de- 
cuplet contributions is crucial to mainlining the correct 
Nc counting as well as the convergence properties of the 
chiral expansion through 0(p2). 

In a recent paper [|3|, we have calculated the 0{p^) 
corrections to J^g arising from octet baryon intermedi- 
ate states. These corrections are entirely of recoil order, 
scaling as inverse powers of the baryon mass. In that 
study, we employed baryon chiral perturbation theory 
with infrared regularization [Q , which effectively resums 
an infinite tower of recoil corrections. Although this re- 
summation is necessary to maintain the analytic proper- 
ties of the currents for momenta near physical thresholds, 
we found that for = the sum is dominated by the 
leading correction which can be obtained directly 

in HBCPT. We also found that the 0{p^) corrections 
were large, exacerbating the poor convergence obtained 
through 0{p^) in octet-only calculations. We left open 
the question as to the impact of including decuplet in- 



termediate states, speculating that large- iVc symmetries 
can generate cancellations at this order as well. 

In the present note, we report on an explicit calcula- 
tion of the 0{p^) corrections which includes contributions 
from the decuplet. We find that - even under the sym- 
metry constraints imposed by the large- A'^c expansion - 
these corrections are both substantial and devoid of the 
cancellations arising at 0{p^). In several channels, the 
0{p^) corrections can be as large as the 0{p'^) term, in 
contrast to the naively expected power suppression by 
(to_r-/A^)2 X {ttik/M) 1/8. We also show that the 
reduced order in Nc arising from the 0{p^) spin-flavor 
algebra is, in retrospect, what one might expect from 
the Nc behavior of the relevant counterterms. In con- 
trast, the 0{p^) loop corrections are finite and entirely 
non-analytic (in quark mass), so there exists no coun- 
terterm at this order whose Nc behavior would imply 
a corresponding order in Nc for the 0{p^) loop correc- 
tions. While this observation does not by itself explain 
the apparent breakdown of the chiral expansion for 
at 0{p^), it does suggest that inclusion of decuplet in- 
termediate states is not generally sufficient to maintain 
the proper scaling behavior of the expansion. As a prac- 
tical corollary, we also note that the use of SU(3) chiral 
perturbation theory to extract As - the strange quark 
contribution to the nucleon spin - from polarized deep 
inelastic scattering data is subject to uncontrolled ap- 
proximations and, therefore, untrustworthy. 



II. AXIAL CURRENTS 

In writing down the octet axial currents, it is con- 
venient to start with the relativistic meson-baryon La- 
grangian. At the lowest order, one has 

Co = iTt {Bijt'D^ - 771n)B) +DTt {Bj^jr^{A^, B}) 
+ FTt (B7'^75[A^, B]) + ift'YDuT^ - uitT^T^ 
+ C[f^'A^B + BA^,T^'] + nT^'Yl5A''T^ 
+ :^Tr((Z?^E)ti?^E)+aTrAf(S + Et)^ (1) 
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One may obtain vector and axial vector current opera- 
tors from Cq by including vector and axial vector sources 
in the covariant derivatives. The leading [C'(p'^)] operator 
contains only baryon fields and the SU(3) reduced matrix 
elements D and F. Axial currents involving both baryons 
and mesons first appear at 0{p). Additional purely bary- 
onic axial currents appear at 0{p^) ||f^. They arise from 
the SU(3) symmetry breaking (SB) Lagrangian 

Ci = ^{di Tr(B7A'75{A^,x+}S) 
+ d2 Tr (B7^75A,,Bx+) 
+ d3 Tr (B7*'75X+BA^) 
+ d4 Tr(B7'^75B{A^,X+})}, (2) 

where 

X = 

Using £0,1 one obtains the axial current: 

= \D Tr (S7^75{^r^C^ + ^T^, B}) 
+ \F Tr (57^75 KT^ft + ^ty^^, b]) 

+ Tr (B7^75{CT^^t + ^ty^^, ^^}b) 

+ id2^ Tr (57^75 (^T^e^ + S}T^i)Bx+) 
+ idg^ Tr {B^^^i^x+BiiT^e + ^^^^0) 

+ id4^ Tr (B7'^75i?{eT-4e^ + ^tT^^, ;^+}) 




i Tr(B7^[eT^et-^tT^e,S]) 



(3) 



The heavy baryon expansion of £0,1 ^nd is ob- 
tained by defining the heavy baryon field H(x) = 
exp{imNV ■ x)^^Y^B{x) [v^ is the baryon velocity) and 
projecting out the postive energy states as in |3[|. In 
this case, all baryon mass terms are removed from the 
Lagrangian at leading order, leaving only a dependence 
on the octet-decuplet mass splitting, 5 = toa — tojv. At 
subleading orders, there are recoil corrections in the form 
of l/ruN- In order to consistently include the decuplet 
we follow the small scale expansion proposed in ||l^]. In 
this approach the energy and momenta and the decuplet 
and octet mass difference 5 are both treated as small ex- 
pansion parameters in chiral counting. Note also that in 
the heavy baryon expansion, one makes the replacement 
7/i I'/j, 7/j75 25"^ etc, where S*^ is the spin operator. 

Renormalized matrix elements of J^^ between octet 
baryon states up to 0{p^) may be written as 



X = TT,K,7J 



+ 



(4) 



where the first term on the right hand side is the lowest 
order one. The second term arises from the SB terms 
in Eq. (||). The third term in Eq. (Q) arises from the 
wave function renormalization. The fourth term comes 
from the vertex correction diagram. The fifth term is 
the vertex correction from the tadpole diagram. The last 
term in Eq. (|^) arises from the 0{p) one- meson operators 
in J^. Details of the last three terms can be found in . 

Terms with Q;^,a^,A^, (if^, and 9^ arise from 
the contribution of octet states only and their expres- 
sions are given in [p^ . The remaining terms come from 
the insertion of decuplet states in the loop. The exprc& 
sions of A:^ '^^ '^^ 
respectively. 

The functions etc are defined as 



■ij > Pij > Pij are presented in Tables |, |l[ and |ll 



— — ) ln( ) -fTT 



*Here, ^ ~ 1 GeV denotes the renormalization scale. In our 
previous analysis ^|, this scale was effectively set equal to 
mjv. Moreover, in that work, the variable ^ denoted the ratio 
niTr/mN and not the renormalization scale. 
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For X — K,ri we have 
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where q and are SU(6) quark creation and annihilation 
operators and A° and cr' are the Geh-Mann and Pauh 
matrices, respectively. At leading order in 1/Nc, one has 



TO 

•^i5 



A" cx G*' 



(12) 



where the coefficient of proportionality is of order unity 
and where terms of relative order 1/Nc have been 
dropped. The Nc counting rules give G" ~ Nc- Thus, 
the ©(p") current is 0{Nc), while loop corrections, which 
contain three insertions of A''"- divided by ^ Nc are 
nominally of order N^. However, the SU(6) commutator 
algebra 



(13) 
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in Eqs. (M) by 



In 



5 + 



TOX 



(8) 



we obtain expressions for f g- I^i t^iis work we explic- 
itly keep pion loop contribution. If we truncate at or- 
der 0(p2) and ignore the pion loops and take S ~ 
and to: 



|to2 



K' 



we reproduce the expressions in 
exactly. Note that we retain only loop corrections hav- 
ing non-analytic dependence on quark masses. Analytic 
terms (e.g, cx to^) have been absorbed into the countert- 
erms di-4. 



implies that the 0{p^) loop corrections, which depend 
on double commutators of are actually of order N^, 
since each commutator reduces the naive counting by one 
power of Nc- 

Because the 0{p^) loops are divergent, there must ex- 
ist counterterms of the same order which absorb the in- 
finities. The most general 0{p^) operators arising at 
this order include those proportional to c?i...4 in Eq. (|^) 
0. These operators involve one insertion of times 
TOp/A^, where mp is the Goldstone boson mass. The 
latter is 0{N°) whereas A^ = {AirF^'f is 0{Nc). Thus, 
the 0(p2) counterterms are 0{N^). Self-consistency of 
the theory implies that the 0{p^) loop corrections must 
also be of 0{N^). Otherwise, there would exist a mis- 
match between the divergent loops and the counterterms 
which render them finite in the large limit. In retro- 
spect, then, one might have anticipated the existence of 
a large- TVc spin-flavor algebra whose affect is to reduce 
the nominal Nc order of the loops to match that of the 
counterterms. 

In contrast, there exist no counterterm operators at 
0{p^), and the loop contributions of this order are en- 
tirely finite and non-analytic in ruq. Thus, one has no 
self-consistency requirement at 0{p'^) involving countert- 
erms and loops to force a reduction in the nominal Nc 



III. Nc COUNTING 

As discussed in a beautiful series of papers |^-[l^] , the 
baryon axial currents have an expansion in 1/Nc involv- 
ing SU(6) spin-flavor operators: 



^ There exist additional operators proportional to and 
as well. The finite parts of the former are numerically 
insignificant while the latter do not contribute to the = 
currents. Thus, we do not show them explicitly, though their 
presence is required to remove the divergences. 
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order of the latter. In particular, the 0{p^) wave function 
renormalization and vertex corrections involve three in- 
sertions of A''"' divided hy x tun. Since niN is 0{Nc), 
these loop effects are nominally order Nc- In the ab- 
sence of any algebra which reduces this nominal order, 
one might expect them to be numerically significant. As 
a practical matter, we find that inclusion of decuplet in- 
termediate states produces no cancellations indicative of 
an algebraic reduction in the nominal Nc order of these 
graphs. Similarly, the 0{p^) seagull graphs involving 
the chiral connection -which have nominal chiral order 
C'(A'^)-receive only octet contributions, so no cancella- 
tion are possible in this case. We also find these con- 
tributions can be significant. Indeed, as we show below, 
the 0{p^) contributions are generally as large or larger 
than the 0{p^) terms, in accordance with naive scaling 
arguments. 



IV. NUMERICAL ANALYSIS 

In tables IV- VII we present various fits to the octet ax- 
ial currents, showing the contributions arising at various 
orders in p. For notational simplicity, we define the axial 



couplings as 



(14) 



where we have omitted the induced pseudoscalar terms. 
In general, we have eight low-energy constants (LEC's) 
to be determined: D, F, (ii_4, TL and C. However, there 
exist experimental data for only six octet matrix elements 
[^5|. Consequently, we must invoke additional assump- 
tions in order to complete the analysis. 

The constants C and 7i can be treated using one of sev- 
eral approaches. Drawing entirely on experimental data, 
the magnitude of C can be determined from the decay 
width of the A. At leading order, one has \C\ = 1.5 [16|, 
which is consistent with the large prediction pO|, 1 1 1 . 
Loop corrections to this result arise at 0{p^). Since C en- 
ters the axial currents at 0{p^), chiral corrections to the 
value of C as determined from the A decay width affect 
our analysis at 0{p'^). Unfortunately, the phase of C can- 
not be determined in this manner, and so one must rely 
on auxiliary considerations. For example, SU(6) symme- 
try implies C and D have the opposite phase. In what 
follows, we make this choice for the phase. 

The situation regarding Ti. is more problematic. This 
LEG does not appear at leading order in any physi- 
cal decay amplitude. It does, however, give the strong 
ttAA coupling at leading order ||l^. A determination 
of this constant is, therefore, highly-dependent on model 
assumptions. In the large Nc limit, for example, Ti. — 
— | (D -t- F). Various quark models yield the same result 
|jl8|-p0|. On the other hand, a light cone QCD sum rule 
analysis |jl^ yields \H\ — 1.35, which is only of half of 
large Nc or quark model prediction and is approximately 



the same value as extracted from from the isobar pro- 
duction experiments in 7r~p iT^'!T~n near threshold 
This constant has also been extracted from decay 
widths using HBCPT to ©(p^) |||. Recently, H was de- 
termined from a fit to phase shift data in the fourth order 
chiral perturbation theory analysis |p3|] . The results im- 
ply 0.94 <n < 2.65. While the magnitude of TL for this 
range is consistent with both the large Nc and QCD sum 
rule analyses, the phase differs from all other approaches. 
It was emphasized in Ref. |23 , however, that H enters 



pion nucleon scattering at third order loop so it can't be 
pinned down precisely. Fortunately, in the case of the 
axial currents, H arises at 0{p^), so the impact of un- 
certainty in this constant is not as pronounced as in the 
case of C. 

A final possibility for treating C and H is to follow the 
analysis of Ref. |q-[l^ and invoke the SU(6) relations: 
C = —2D, Ti = — Doing so reduces the number of 
fit parameters to six |^ The authors of Refs. |]8|-|l2| found 
that use of SU(6) relations among the LEC's minimizes 
the size of the 0{p'^) loop corrections, in accordance with 
the cancellations expected from large Nc arguments. It 
is not possible to apply similar relations to di-4, how- 
ever, since they parameterize explicit symmetry-breaking 
terms in the Lagrangian. 

In Table IV, we give a fit through 0{p^) using the 
experimentally-determined magnitude for C, a phase op- 
posite to that of D, and the quark model value for H. 
The remaining six LEC's are determined from the nu- 
cleon and hyperon semileptonic decay data. Under these 
conditions, the 0{p^) corrections (including both loop 
effects and symmetry breaking terms) are generally as 
large as the 0{p^) contributions. However, invoking the 
SU(6) relations among D, C, and H changes this situa- 
tion considerably, as illustrated in Table V. In this case, 
the relative importance of the 0{p^) terms is consider- 
ably reduced and the improved. In Table VI, we show 
the corresponding fit using the SU(6) relations but set- 
ting di_4 — 0. The latter fit corresponds roughly to the 
analysis of Ref. |p|-p^, which illustrated the impact of 
O(p^) loop cancellations in the symmetry limit. Gener- 
ally speaking, inclusion of di-4 improves the quality of 
the fit as well as the scaling behavior of the chiral expan- 
sion through O(p^). 

In Table VII we give the best fit through 0{p^). Here, 
we have used the SU(6) relations for D, C, and H in order 
to produce the cancellations at 0{p'^). We observe that 
the 0{p^) contributions are generally as large or larger 



"*We have also used the relations arising from the inclusion 
of 1/iVe corrections to Eq. (|l|): C = -2D, H = 3D - 9F in 
our fit. The fit results turn out to be the same. We thank E. 
Jenkins for suggesting this point. 

further reduction in number of parameters may occur 
when the double expansion in and l/Nc of Ref. |Q is 
used. 



4 



than the C(p^) terms and, in several channels, as large 
as the 0{p^) terms. This pattern becomes even more 
pronounced away from the SU(6) limit for the LEC's, in 
which case neither the 0{p^) nor the 0{p^) terms scale 
as expected. 

The breakdown of the chiral expansion which wc ob- 
serve at at 0{p^) reflects a number of factors: the large 
magnitude of the kaon mass, which appears in the numer- 
ator of the expressions in Eqs. (||); the apparent absence 
of cancellations (and an underlying large Nc spin-flavor 
algebra) among the recoil order corrections; and the ap- 
pearance of factors of tt in integrals and arising 
at this order. 



V. DISCUSSION 

It has been known for many years that tree-level SU(3) 
relations are remarkably successful in describing a num- 
ber of the low-lying properties of hadrons, such as pseu- 
doscalar masses and baryon axial currents. Ideally, chiral 
perturbation theory - together with the large expan- 
sion ~ should suffice explain why these relations work 
so well. With such an understanding in hand, one would 
have had considerable confidence in exploiting these rela- 
tions to determine quantities for which one has no direct 
measurement, such as the strange quark contribution to 
the nucleon spin. As. In the present study, however, we 
observe that the chiral expansion for baryon octet axial 
currents does not appear to be under control. While large 
Nc considerations imply that the expansion works rea- 
sonably well through ©(p^), it breaks down completely 
at 0(p^)p*|. While a theoretical justification for applying 
SU(3) symmetry to the octet axial currents may exisij^, 
we are unable to provide one at this time. 

As a practical consequence of this situation, we con- 
sider the determination of As from polarized deep inelas- 
tic scattering (DIS) data. As shown in Ref. [Q, one may 
express As in terms of the polarized structure function 
integrals 



1 



dx 5i'"(a;) 
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The numerical separation of through 0{p^) is given 
in Table VII and yields 



As = 0.14 - [0.12 + 0.25 -I- 0.10] 



(17) 



where the numbers in square brackets correspond, re- 
spectively, to the order p", p^, and p^ contributions to 
gg. Since the chiral expansion is not converging for As, 
we do not quote a total for this quantity nor can we 
estimate a theoretical uncertainty. In contrast, extrac- 
tions of As from semi-inclusive measurements performed 
by the Hermes collaboration ||2^ or from elastic neutrino- 
nucleon scattering ||2^,^ are not plagued by large SU(3)- 
breaking uncertainties, making them in principle more 
reliable probes of the flavor content of the nucleon spin. 

This work was supported in part under U.S. Depart- 
ment of Energy contract #DE-FGO2-00ER41146, the 
National Science Foundation under award PIIYOO-71856, 
and a National Science Foundation Young Investigator 
Award. We thank E. Jenkins, B. R. Holstein, and J. L. 
Goity for useful discussions. 



As = -[rp 



r„j - —98 



(16) 



where g^ is the axial vector coupling associated with the 
matrix element {p\J^\p). The combinations of LEC's re- 
quired for this matrix element are 



** However, we observe that at C'(p^), there exist some chan- 
nels for which the large Nc cancellations are not strong (see, 
e.g., Table 0). 

^^See, e.g., the regulator scheme proposed in Ref. [M 
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TABLE II. The coefficients (3^^ for the vertex corrections 
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TABLE 111. The coefficients for the vertex corrections. 
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TABLE IV. The separation of fit results into pure 
©(p") and 0{jP) pieces where we have used 5 = 0.3 GeV, 
mjv ^ 00, C = —1.5, Ti = —2.25 as inputs. The fit yields 
D = 0.63, F = -0.45, di = 0.79, d2 = 1.87, = 1.43, d4 = -1.13 
with = 0.15. 
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Full fit results 


Tree level only 


0(p'') only 


.jpn 


i.2() 


0.77 


O.i!) 


A 
Aa, 


0.62 


0.38 


0.24 


9pA 


-0.89 


-0.57 


-0.32 


5ab- 


0.32 


0.19 


0.13 




0.34 


0.15 


0.19 


5^0 5- 


0.92 


0.54 


0.38 




0.15 


0.15 







0.17 


0.13 


0.04 



TABLE V. The separation of fit results into pure 0{p°) 
and 0{p^) pieces where we have used 5 = 0.3 GcV, 
mjv — > oo, C = —2D, TC = —3D as inputs. The fit yields 
D = 0.46, F = 0.31, di = -0.80, da = 0.93, ds = -0.63, = 0.78 
with = 0.002. 





E\ill fit results 


Tree level only 


0{p') only 




1.10 


0.76 


0.34 




0.66 


0.42 


0.24 


9pA 


-0.88 


-0.51 


-0.37 




0.31 


0.10 


0.21 


A 
9nS- 


0.29 


0.26 


0.03 


ffEOS- 


1.05 


0.54 


0.51 


9io^-' 


0.35 


0.26 


0.09 




0.26 


0.07 


0.19 



TABLE VL The separation of fit results into pure 
and 0{p^) pieces where we have used 5 = 0.3 GeV, mjv — > oo, 
di-4 = 0, C = -2D, n = -3D as inputs. The fit yields 
D = 0.51, F = 0.25 with = l-l- 





Full fit results 


Tree level only 


0{p') only 


0{p') only 


9pn 


1.26 


0.61 


0.41 


0.24 


3ae- 


0.58 


0.32 


0.14 


0.12 


9pA 


-0.92 


-0.43 


-0.11 


-0.38 


9L 


0.26 


0.11 


0.05 


0.10 




0.33 


0.17 


0.03 


0.13 


5eos- 


0.87 


0.43 


0.05 


0.39 




0.22 


0.17 


-0.02 


0.07 




0..32 


O.OS 


0.17 


0.07 



TABLE VIL The separation of fit results into pure 0{p°), 
0{p^), and 0{p^) pieces where we have used 5 = 0.3 GeV, 
ruN = 0.94 GeV, C = -2D, H = -3F as inputs. The fit 
yields 

D = 0.39, F = 0.22, di = -1.97, di = 1.14, ds = -0.45, d^ = -0.06 
with x^ = 0.12. 
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